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ABSTRACT 



Here we attempt to infer the recent history of star formation in the BCD galaxy VII Zw403, 
based on an analysis that accounts for the dynamics of the remnant generated either by an 
instantaneous burst or by a continuous star formation event. The models are restricted by the 
size of the diffuse X-ray emitting region, the luminosity from the star-forming region and the 
superbubble diffuse X-ray luminosity. 

We have re-observed VII Zw403 with a better sensitivity corresponding to the threshold Hq, 
flux 8.15 X 10~^^ erg cm~^ s~^. The total Hq, luminosity derived from our data is much larger than 
reported before, and presents a variety of ionized filaments and incomplete shells superimposed 
on the diffuse Hq, emission. This result has a profound impact on the predicted properties of 
the starburst blown superbubble. Numerical calculations based on the HST Hq data, predict 
two different scenarios of star formation able to match simultaneously all observed parameters. 
These are an instantaneous burst of star formation with a total mass of 5 x 10'"' M© and a star- 
forming event with a constant SFR = 4 x 10"'^ Mq yr"-"^, which lasts for 35 Myrs. The numerical 
calculations based on the energy input rate derived from our observations predict a short episode 
of star formation lasting less than 10 Myrs with a total star cluster mass (1 — 3) x 10^ Mq. 
However, the five main star-forming knots are sufficiently distant to form a coherent shell in a 
short time scale, and still keep their energies blocked within local, spatially separated bubbles. 
The X-ray luminosities of these is here shown to be consistent with the ROSAT PSPC diffuse 
X-ray emission. 



Subject headings: galaxies: starburst 
bubbles — ISM: abundances 

Introduction 



galaxies: dwarf — galaxies: individual (VII Zw403) — ISM: 



It has recently been recognized that the star 
formation activity in galaxies is very irregular in 
time, and many examples of major burst episodes 
exhibit an extremely high star formation rate con- 
centrated in well localized space regions (Terlevich 
1996). It is also now well known that starbursts 
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(SBs) cause an emission that dominates the en- 
tire host galaxy luminosity and their mechani- 
cal energy input rate is expected to cause major 
structural changes in the surrounding interstellar 
medium (ISM). In this respect it has become of 
great interest to study the properties of the re- 
sultant large-scale expanding superbubbles which, 
powered by the violently injected newly processed 
matter, establish the time scale for mixing with 
the ISM (Tenorio-Tagle 1996, Sifich et al. 2001). 
In extreme cases, the superbubbles are thought to 
break out of the galactic discs leading to an effec- 
tive mass and energy transport into the low den- 
sity halos or even into the intergalactic medium 
via a superwind (Heckman et al. 1990). 

Starbursts in the local universe are also as- 
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sumcd to be good representatives of the star- 
forming activity at high redshifts. This con- 
cept defines their cosmological interest as key 
laboratories for studying the ISM, the transport 
of supernovae processed metals, as well as the 
chemical evolution of galaxies and of the inter- 
galactic medium. The resulting structure in the 
ISM due to mechanical energy injected by SBs is 
very similar to the interstellar wind-blown bub- 
bles around single massive stars (sec Weaver et 
al. 1977 for their four zone model), although the 
much larger energy input rate in SBs leads rapidly 
to much larger scales. Hydrodynamical simu- 
lations (see Tenorio-Tagle & Bodenheimer 1988, 
Bisnovatyi-Kogan & Silich 1995 and references 
therein; Suchkov et al. 1994; Silich & Tenorio- 
Tagle 1998; D'Ercole & Brighenti 1999; Strick- 
land & Stevens 2000) currently include differential 
galactic rotation, radiative cooling, strong density 
gradients between the disk and the halo and thus 
are able to follow the moment of breakout, as well 
as the fragmentation of the expanding outer shell 
via Rayleigh - Taylor instabilities and the venting 
of the super bubble hot interior gas, either into the 
intergalactic space, or into the host galaxy halo 
(the blowout phenomenon). 

Most of the up to-date simulations have been 
performed under the assumption of a constant en- 
ergy deposition rate, as expected from an instan- 
taneous burst model. However, studies of the stel- 
lar population in OB associations related to young 
(roB < 10 Myr) Large Magellanic Cloud (LMC) 
bubbles (see Oey & Smedley, 1998 and references 
therein) have demonstrated that "realistic" energy 
input rates are very different from the assumed 
constant energy input rates used in numerical sim- 
ulations. Thus the instantaneous burst assump- 
tions may not be applicable to all cases. Here we 
attempt to establish a method of comparison be- 
tween the theory of superbubbles and the observa- 
tions of remnants produced by massive star forma- 
tion in galaxies. Two possible modes of star forma- 
tion, instantaneous and extended bursts, are taken 
into consideration. For both cases the mechan- 
ical luminosity, ultraviolet photon output, mass 
returned to the ISM, and the fraction of each in 
metals, all as a function of time, are estimated. On 
the other hand we have the observed parameters: 
the Hq, or Hp luminosity which can be directly 
related to the SFR under the assumption that all 



photons are used up in the ionized region. One can 
also estimate the size and luminosity of the X-ray 
remnant. In some cases the remnants may have 
slowed down sufficiently to display their outer ex- 
panding shells, either in the optical or in HI obser- 
vations, or both, giving further information about 
the size, expansion speed, and mass behind the 
outer shock. A comparison with the theory also 
requires some preconception of the galaxy's ISM, 
which can as a first approximation be derived from 
HI observations and the inferred dynamical mass. 
However, one also needs to make an assumption 
about the fraction of this gas locked up in dense 
clouds and immersed into a less dense ISM back- 
ground. 

With the aim of establishing a method to con- 
front theory with observations that may lead us 
to infer the recent star formation history of galax- 
ies, the blue compact dwarf (BCD) VII Zw403 is 
thoroughly analyzed. Section 2 derives the main 
properties of coeval and extended star formation 
modes. Section 3 summarizes the main observa- 
tional properties of our target. Section 4 presents 
a summary of the hydrodynamical calculations 
aimed at matching the observed properties. The 
results of the calculations and our main findings 
are discussed in Section 5. 

2. Star formation 

Following the observational evidence in sup- 
port of the concept of a continuous star formation 
spread over 10 - 20 Myr or more (Herbst & Miller 
1982, Stabler 1985, Coziol et al. 2001), Shull & 
Saken (1995) have discussed five possible tracks of 
the star formation rate (SFR) time evolution, in- 
cluding both the time span of star formation, and 
possible time variations in the initial mass func- 
tion (IMF). Here we examine one of their cases: 
a constant SFR spread over a finite time interval 
TSB- We compare its intrinsic properties with the 
standard instantaneous model. The simple model 
discussed below can be generalized to more sophis- 
ticated cases with a time dependent SFR and a 
variable IMF. 
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2.1. The properties of starbursts 

Let us assiimc an instantaneous burst charac- 
terized by an IMF with an exponent a, 
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where Msb is the total mass of the SB cluster, 
and Miow and M„p are the lower and upper cut-off 
masses, respectively. From this one can calculate 
the expected energy and mass input rates. 

2.1.1. The mechanical energy input rate. 

The mechanical energy deposition from a star- 
burst should include both the energy injection via 
stellar winds and from supernova explosions. Here 
the early mechanical energy injection (before the 
first supernova explosion a,t t = tgn) has been ap- 
proximated by a constant value L^, which adds 
the contribution of all individual massive stars in 
the cluster. This value can easily be scaled from 
the template of Leitherer & Heckman, 1995 (here- 
after LH95) for a starburst with a total mass of 
10^ Mq within the mass range Miow = 1 M© and 
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where Msb is the total stellar mass considered, 
Mlh — 10® M© and L^h, following LH95 was set 
equal to 10^^ erg s~^. 

After a time {t > tsn) the energy released via 
supernova explosions dominates the SB mechani- 
cal energy deposition, and then the energy released 
by the SB as a function of time is defined by the 
expression 

EsB{t)= / EsNf{m)dm = 

JM{t) 



,(3) 



where M{t) is the mass of the stars exploding as 
supernovae after an evolutionary time t. For sim- 
plicity, the energy released per supernova {Esn) is 
assumed to be independent of the progenitor mass 
and equal to 10^^ erg. 



The energy input rate from the SB at the su- 
pernova dominated stage is then 

dEsB{t) 
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The lifetime of the massive stars can be inferred 
from the approximations of Chiosi et al. (1978) 
and Stothers (1972). Then the function M{t) has 
the form: 

M{t) J (9 X 10Vt)2; 3OM0 < M < IOOM0 
IOMq ~ [ (3 X 107t)5/8; 7Mq <M < SOMq. 

(5) 

Substituting the time derivative of M{t) into equa- 
tion (4) one obtains 
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where tso and tr are the stellar lifetimes for a 30 
M and a 7 M © star, respectively. 

2.1.2. Mass deposition 

Before the first supernova explosion, the mass 
injected by a star cluster results from individual 
stellar winds calculated as 

Mej(i) = / 2LsB/V^dt, t < tsn, (7) 



where the collective wind terminal velocity is as- 



sumed to be V,,, = 2000 km t 
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Afterwards, 



the mass injection is dominated by supernovae 
and one can neglect the contribution from stel- 
lar winds. During the supernova dominated stage, 
the total ejected mass as a function of time may 
be approximated by (see Silich et al. 2001): 
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Following Silich et al. (2001) we assume that 

the gas ejected by supernovae includes all the 
newly synthesized metals, and thus the mass in 
ejected metals as a function of time is 
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where YmetaU is a particular clement yield. One 
can then use, for example, oxygen as tracer of the 
metalicity caused in the interior of super bubbles. 
In such a case the oxygen yield can be approxi- 
mated by analytical fits to the stellar evolutionary 
tracks of Maeder (1992) and Woosley et al. (1993), 
which account for the mass loss due to stellar 
winds (see Silich et al. 2001). The mean metalicity 
inside a superbubble is then defined by the ratio of 
the wind and supernovae ejected oxygen to the to- 
tal mass within the superbubble Mj„ = Mej+M^v, 
which includes also the mass thermally evaporated 
from a cold outer shell {Mev)- 

M^j{0)/ZQ{0) + ZisMM^^ 



Zo 



Mi, 



(10) 



where Zjsm and Zq are both in the solar units, 
Zq{0) = 0.0083 (Grevesse et al. 1996). 

2.1.3. The UV output 

For a coeval star cluster, the UV photon pro- 
duction rate remains almost constant until the 
most massive star in the cluster begins to move 
away from the main sequence, soon becoming a 
supernova. From this moment onward, the UV 
flux rapidly decays (Beltrametti et al. 1982). We 
approximate this evolutionary track by the power 
function 



Nuv 
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The initial flux {N^) was normalized to the 10^ 

M0 standard model of LH95, and was taken to 
be 9 xlO^^ photons s~^. If all UV photons are 
trapped within a star forming region, then the 
H„ expected luminosity L}j^ follows from a simple 
transformation (LH95): Lh^ = 1.36 x IQ'^'^Nuv 
erg s~^. 

2.2. Continuous star formation 

If the star formation process is not coeval but is 
instead spread over a time tsb, one can approxi- 
mate it by a series of Ntot instantaneous sequential 
mini bursts separated by a At = Tss/Ntot- One 
can further assume that all mini bursts would have 
the same mass Mi = Msb /Ntot and IMFs and 
would evolve independently according to their own 
clocks set upon formation: 



where t is an evolutionary time, and i is the mini 
burst number. Then the number of mini bursts at 
any given time t is 

r Int{t/AT) + 1, if t<TSB ^^^^ 

n{t) = < (13) 

t Ntot, tj t> TSB- 

The cumulative intrinsic properties (mechanical 
energy input rate, mass injection and the number 
of UV photons or the Ha luminosity) can be found 
by simply adding the mini burst parameters: 
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Note that the input from each mini burst should 

be set to zero whenever their individual evolution- 
ary time ti exceeds the lifetime of the less massive 
star which can explode as supernova; that is once 

ti becomes larger than their intrinsic ty. 

The global energy properties of star forma- 
tion events derived from this simple model are 
in a good agreement with the two extreme cases 
(instantaneous and continuous star formation) of 
LH95 and are presented in Figure 1. For compar- 
ison a value of a = 2.35 has been assumed as well 
as a stellar range of masses between 1 Mq and 100 

Mq. 

Figure la compares the mechanical energy in- 
put rates derived for a 10^ Mq stellar cluster 
formed during a, tsb = 1 Myr (hereafter the in- 
stantaneous burst), 20 Myr and 50 Myr. The solid 
line shows the rapid rise in the coeval case, which 
reaches a maximum immediately after the onset of 
supernova explosions, and then drops to values ~ 
10"^°'^ erg s~^ to remain almost constant up to the 
end of the supernova activity (^ 50 Myr). For the 
TSB = 20 Myr event, the maximum energy depo- 
sition is delayed in time, and converges with the 
1 Myr values after approximately 25 Myr. The 
energy deposition rate for r = 50 Myr star form- 
ing episode suffers an even longer delay, to reach 
values similar to the maximum input from the pre- 
vious two models, and then slowly decays without 
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reaching a uniform value throughout the 100 Myr 
of its evolution. 

The production rate of UV photons (see Figure 
lb) is highly dependent on the star formation time 
considered. For the 1 Myr starburst the UV flux 
reaches its maximum value (~ 9 xlO^^) at t = 
TsB, and remains constant during the first 4 Myr. 
It then falls very steeply, reaching before 10 Myr of 
evolution, values more than two orders of magni- 
tude smaller than its maximum value. The 20 Myr 
and 50 Myr clusters acquire their uniform constant 
value after 4 Myr of evolution, and then upon com- 
pletion of their star formation phase, their UV 
photon output also rapidly decays with time. 

The rate at which mass is ejected by the mas- 
sive star cluster is shown in Figure Ic. The to- 
tal ejected mass 40% of the star cluster mass) 
and the total oxygen mass released by supernovae 
(about 2% of the total mass turned into stars) is 
clearly the same at the end of the evolution in the 
three cases considered. There are however sub- 
stantial delays for the larger the value of tsb- 

3. The main observational properties of 
VII Zw403 

VII Zw403 is a BCD galaxy considered in recent 
years in many discussions related to star forma- 
tion histories and possible impact of dwarf sys- 
tems on the surrounding intergalactic medium. 
Although TuUy et al. (1981) proposed that the 
galaxy might be a member of M81 group, more re- 
cent distance determinations (Lynds et al. 1998) 
locate it about 1.5 Mpc further away at 4.5 Mpc 
distance. Therefore it can be considered an iso- 
lated galaxy (Schulte-Ladbeck et al., 1999). 

The X-ray observations (Papaderos et al. 1994; 
Fourniol 1997) added more interest to the system, 
as they revealed an extended kpc-scale region of 
diffuse X-ray emission. These observations ac- 
count for approximately 85% of the total X-ray 
luminosity {Lx, total ~ 2.3 x 10^* erg s^^) from the 
central unresolved core, while the remaining 15% 
is spread over the central kpc-scale region. The 
high-resolution data reported by Lira et al. (2000) 
reveakxl a strong point X-ray source displaced to 
the west with respect to the main Hq, emission. 
This is most probably a powerful (Lx ~ 10^^ erg 
s~-^) binary system. However, the kpc-scale dif- 
fuse component was not confirmed by these obser- 



vations, possibly due to the low sensitivity of the 
HRI. 

The total dynamical mass of the galaxy is 2 x 
10^ Mq, with approximately 20% or 4 x lO'^ M© 
constituting the neutral hydrogen phase (Thuan 
fc Martin, 1981). VII Zw403 does not exhibit a 
regular rotation, but rather turbulent random mo- 
tions of individual neutral hydrogen clumps with 
AVt w 30 km s'^ (Thuan & Martin, 1981). This 
value is taken as the velocity dispersion which 
supports the ISM against the self gravity of the 
galaxy. The extension of the neutral hydrogen 
halo RjsM is not well known, but if the average 
ratio of the HI size to the optical size is accepted 
to be f = 2.4 (Thuan & Martin, 1981), the HI halo 
radius Rhi is about 1.9 kpc. A similar result is 
obtained if one assumes that the galaxy boundary 
occurs at the location where the escape velocity 
drops below the ISM random velocity AV*. 

Lynds et al. (1998) showed that the stellar pop- 
ulation of VII Zw403 contains not only young blue 
main-sequence stars, but also a much older evolved 
population that can be traced up to 1-2 Gyrs back 
in time. They also studied the five centrally con- 
centrated associations of young stars, with ages 
smaller than 10 Myr, and that in Hq, emission ap- 
pear to be surrounded by local small shells with 
radii ~ 100 pc. Schulte-Ladbeck et al. (1999) pro- 
vided a near infrared single star photometry of the 
galaxy and confirmed that compact, young star- 
forming regions are embedded into a much older, 
low surface brightness halo. 

3.1. New observations 

We present here new narrow-band images of 
VII Zw403 centered on the Ha line and on the ad- 
jacent continuum (Figure 2) which were taken 
in August 1997 at the 2.2-m telescope of the 
German-Spanish Astronomical Observatory at 
Calar Alto (Almerfa, Spain). The instrumen- 
tation consisted of the Calar Alto Faint Object 
Spectrograph (CAFOS) and a 2048 x 2048 SiTe 
GCD chip, with a pixel size of 0.53 arcsec, and 
an un- vignetted circular field of view of about 11 
arcmin in diameter (14.4 kpc at the accepted dis- 
tance to the galaxy of 4.5 Mpc). The averaged 
seeing was 1.5 arcsec. 

The image reduction was conducted using stan- 
dard procedures available in IRAF. Each image 
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was corrected for bias using an average bias frame, 
and was flattened by dividing by a mean twilight 
flat field image. After, they were registered (for 
each filter wc took a set of dithered exposures) 
and combined to obtain the final frame, with cos- 
mic rays removed and bad pixels cleaned. The 
average sky level was estimated by computing the 
mean value within various boxes surrounding the 
object, and subtracted out as a constant. Flux cal- 
ibration was done through the observation of spec- 
trophotometric stars from Oke (1990). For more 
details about data reduction and calibration refer 
to Cairos et al. (2001). We calculated the inte- 
grated Hq, flux out to the limiting isophote with a 
level equal to 2.5 x rms of the background. In our 
case this threshold is 8.15 x 10^^''' erg s"-'^ cm~^. 
The Ha flux was corrected for Galactic extinction 
following Burstein & Heiles (1984). No correction 
for internal extinction was performed. 

Our Hq, image shown in Figure 2 displays a va- 
riety of shapes and structures which deserve dis- 
cussion. Besides the starburts knots, labeled 1 - 
5 after Lynds et al. (1998), we detect a much 
fainter and extended well structured diffuse emis- 
sion. The largest obvious feature is a broken to the 
SW shell-like structure, 250 pc in radius, which we 
associate with knot 4. Knot 2 is almost overlap- 
ping in projection along its Eastern rim. 

On top of features and structures which appear 
to build a network of connecting knots, there is 
down to our limiting flux a diffuse-smooth emis- 
sion which engulfs them all and most certainly re- 
sults from photons leaking out of the main emit- 
ting knots. 

The total H„ luminosity is higher than reported 
previously. Clearly our limiting flux is much lower 
than that of Lynds et al. (their Figure 2). Our 
data recovers the low intensity emission and 
thus provides a more precise estimate of the total 
luminosity. 

The spatial resolution of the HST image (Lynds 
et al. , 1998) also allows for the detection of 
a much smaller mini-bubble (cocoon), asso- 
ciated with knot 1. This bubble has similar en- 
ergy to that inferred for the large superbubble re- 
ported above, however it seems much less evolved. 
The energy dumped by the cluster energizing knot 
1 has not been able to sweep and displace so 
efficiently the surrounding ISM out of which it 
formed. 



The above result has a profound impact on the 
modeling of starbursts. Clearly during the flrst 
few Myr of evolution every center of star forma- 
tion disperses its high density parental cloud. It 
is only after this task has been completed that the 
various star forming centers that compose a star- 
burst would be able to jointly build up a large- 
scale superbubble. 

The main observational properties of VII 
Zw403 from the above mentioned literature and 
our own observations are summarized in the Table 
1. 

4. A numerical model for VII Zw403 

A successful model has to be able to match 
simultaneously all observed parameters. In the 
case of VII Zw403 these are: the size [Rsh ~ 1 
kpc) of the diffuse X-ray emitting region, the Ha 
and the X-ray (L^, ~ 3 x 10'^^ erg s^^) luminosi- 
ties. Note that Lynds et al. (1998) give a total 
value IjHa ^ 1.8 x 10^^ erg which is smaller 
than the total ffux detected by our observations 
{Lho. ~ 4.7 X 10^° erg s~^). We have used both 
of these values in our search of a successful model 
of VII Zw403. Note also that the values derived 
from X-ray observations are doubtful and require 
confirmation (Bomans, 2001). VII Zw403 presents 
various nuclear centers of star formation, each of 
which have managed to locally structure the ISM, 
and thus present a number of loops or broken shells 
with radii ~ 100 pc. Another important fact is 
that available observations display neither an HI 
nor an HII or X-ray large-scale shell surround- 
ing the diffuse X-ray emission. This fact can be 
ascribed to a remnant evolving along its quasi- 
adiabatic phase, or to a neutral hydrogen shell 
brightness that falls below the detection limit. 

All considered models assume that the observed 
HI mass occupies a smooth low density disc-halo 
distribution, although an important fraction of it 
(20% < fc < 95%) is in a dense cloud compo- 
nent. This has a major impact on the evolution 
of remnants as it affects both the time required to 
reach a given size, as well as the resultant X-ray 
luminosity produced within the superbubble. Our 
models also assume that the ISM is supported in 
hydrostatic equilibrium in the galaxy gravitational 
field by the turbulent gas pressure produced by 
the random gas motions with an effective velocity 
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dispersion AVt (sec Tomisaka & Brcgman, 1993; 
Silich & Tenorio-Tagle, 1998). An example of the 
resultant initial ISM smooth component density 
distribution is shown in Figure 3. 

The calculations were carried out with our 3D 
Lagrangian code, which accounts for the enrich- 
ment of the hot superbubble interior by the met- 
als ejected via supernova explosions (Silich et al. 
2001). Using oxygen as tracer, we derived the 
time dependent superbubble interior gas metal- 
icity, which was then used in our calculations of 
the diffuse X-ray emission. As there are no firm 
observational restrictions on the VII Zw403 star 
formation timescale, we have examined different 
scenarios of the recent star formation activity in 
this galaxy. Two sets of models have been con- 
sidered: an instantaneous burst, and an extended 
star formation episode. We have associated an 
instantaneous burst with the 1 Myr event and cal- 
culated the SB parameters following the prescrip- 
tions of section 2. In all coeval models consid- 
ered, the energy input rate has been derived from 
the assumed mass of the star forming cluster. In 
the extended star formation scenarios we used in- 
stead the observed H„ luminosity to derived the 
star formation rate from our continuous star for- 
mation model (see section 2). This SFR was then 
transformed into a total star cluster mass and en- 
ergy input rate, for the various assumed values of 
TSB, as it was considered in section 2. A Salpeter 
initial mass function, with upper and lower cutoff 
masses 1 and 100 M© respectively, was assumed 
for all models. 

We then calculated the time Tdyn that super- 
bubbles need to reach the observed 1 kpc diffuse 
X -ray radius. A Vciluc of Tdy^i implies an age from 
which one can calculate the central H„ luminos- 
ity and the integrated superbubble diffuse X-ray 
emission for every model. The results of the cal- 
culations for instantaneous burst models are sum- 
marized in Table 2. The various models are la- 
beled with two separate indices that represent the 
logarithm of the SB mass (also indicated in col- 
umn 2) and the fraction of the ISM assumed to 
be in the smooth component respectively. The 
preceding letter "1" indicates that all of these are 
instantaneous burst models. Four different star 
cluster masses have been considered: 10^, 5 x 10^, 
10® and \{fM.Q. Column 3 indicates the fraction 
of the ISM mass assumed to be stored in clouds. 



The calculated dynamical time Tdyn is presented 
in column 4. After this time the shock wave has 
expanded to Ikpc radius, and the relevant param- 
eters from the model can be compared with the ob- 
served values. The calculated number of UV pho- 
tons, the Ha luminosity (under assumption that 
all UV photons are trapped within a SB region), 
and the resultant superbubble diffuse X-ray lumi- 
nosity, are presented in columns 5, 6 and 7 respec- 
tively. Column 8 indicates whether the superbub- 
bles evolve along a quasi-adiabatic track (A), or 
have made a transition to a radiative phase (R) 
before reaching the 1 kpc radius. The two possi- 
ble modes of expansion imply different properties 
for the swept-up matter: either material is still 
too hot to radiate effectively and thus occupies 
the outer 20% of the remnant volume, or it has 
collapsed into a thin and cold outer shell, while 
being exposed to the ionizing radiation. 

The second set of models assumes an extended 
phase of star formation. The results of the calcu- 
lations for these models are summarized in Table 
3. The various models are labeled again with two 
separate indices that now represent the star for- 
mation time TSB (in Myr) and the fraction of the 
ISM in the smooth component. The letter "C" 
indicates the continuous star forming mode, for 
which episodes lasting 5, 10, 20 and 40 Myr have 
been assumed. In all these models we assumed 
the same SFR (4 x 10"^ M© yr^^). Column 2 
indicates the total stellar mass, and columns 3 to 
8 list the same variables considered in Table 2. In 
all of these models we assumed that the photons 
currently produced by the stellar clusters are all 
used up in reestablishing the ionization of the cen- 
tral HII region. This fact is supported by the large 
HI mass present in VII Zw403. 

The resulting H^ luminosity, and the time de- 
pendent mechanical energy input rate, are shown 
in Figures 4a and b respectively. 

5. Discussion 

Table 2 shows that both low (10^ M©) and high 
(lO'' M©) mass SB models are completely inconsis- 
tent with the VII Zw403 parameters derived from 
Lynds et al. (1998) data. Indeed, the low mass 
models lead to a very slow expansion speed, which 
makes them reach the Ikpc size only after 10 Myr, 
even in the lowest density case with a cloud mass 
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filling factor fc = 90% of the observed HI mass. 
After this time the most massive stars have left 
the main sequence and exploded as supernovae, 
causing a fast drop in the number of emitted UV 
photons and consequently in the luminosity of 
the associated HII region. This luminosity in all 
low SB mass cases is several orders of magnitude 
below the currently observed value. The high mass 
models on the other hand, are too energetic, pro- 
duce an exceedingly high X-ray emission, and an 
overwhelming UV photons flux when the super- 
bubble radius reaches Ikpc. 

The intermediate mass SBs (5 x 10^ - 10^ Mq) 
are in better agreement with observations. How- 
ever, for the emitted number of UV photons to 
be consistent with the value derived from the ob- 
served B.a luminosity, the smooth component of 
the ISM has to contain only a small 5%) frac- 
tion of the observed HI mass. The shock wave 
blown by coherent SN explosions into this low den- 
sity medium remains adiabatic when it reaches a 
1 kpc radius (see Table 2). Therefore, in these 
models (I5.7_5, I6.0_5) an essential fraction of the 
X-ray emission (~ 80% in the model I5.7_5 and 
~ 50% in the model I6.0_5) arises from the outer 
adiabatic shell of swept-up matter that should oc- 
cupy the outer < 20% of the remnant volume. 
However, this is not resolved in the available X- 
ray maps. 

From the results in Table 3, it is clear that 
for a reasonable agreement between the observed 
parameters of VII Zw403 and a continuous star 
formation scenario, the episode of star formation 
should last no less than the dynamical time re- 
quired for the superbubble to reach the 1 kpc ra- 
dius. 

The best continuous star- forming model (C40_70), 
with a SFR = 4 x IQ-^ M0 yr^^, tsb = 40 Myr, 
and 30% of the total ISM mass in the dense cloud 
component, requires about 35 Myr for the outer 
shock to reach a radius ^ 1 kpc. In this model 
an extended gaseous halo cannot completely pre- 
vent gas loss from the galaxy and allows a final 
shell speed slightly in excess of the escape velocity. 
The Hq, luminosity is not an issue in this case as 
it is exactly the amount used to derive the con- 
stant SFR. The time evolution of the superbubble 
diffuse X-ray emission is shown in Figure 5 for 
different ISM models. Initially the X-ray luminos- 
ity grows rapidly with time, following the energy 



input rate and mean hot gas mctalicity time evo- 
lution. It reaches a maximum value after 10 - 20 
Myrs, and then remains almost constant up to the 
end of the star forming activity. The maximum 
X-ray luminosity is larger in models with a larger 
smooth ISM component, and matches the value 
observed in VII Zw403 when the smooth compo- 
nent of the ISM amounts to 70% of the observed 
HI mass. Note that we stopped our calculations 
when the shock fronts reached the assumed galac- 
tic ISM cut-off position (1.9 kpc). This causes 
the breaks in the X-ray curves around 30 Myr 
and 50 Myr, for the fc = 0.9 and fc = 0.6 models 
respectively. 

Note that for these moderate ISM densities and 
the relevant mechanical luminosities (see Figure 
4b), the shell of swept-up matter becomes radia- 
tive after a short time TcooI ~ 2.3 x 10*717^^^^38^^ 
yr (Mac Low & McCray, 1988). However, this 
is not detected in the available HI maps. One 
can then claim that the large-scale shell is photo- 
ionized by radiation escaping the central HII re- 
gion, and thus it should be most easily observed 
in Hq light. However, note that if ~ 80% of the 
UV photon flux is being used to ionize gas around 
the central star clusters, as it is assumed to be 
for the totality of clusters 1 and 2 in VII Zw403 
(Lynds et al., 1998), then the rest of the photons 
(Ne^c = 0.25LHa/l-36 X 10-12 s-i « 3.25 x 10^° 
s^i) would be free to ionize the outer ISM. If the 
number of UV photons trapped within the neutral 
smooth component is negligible, they will ionize 
the neutral outer shell. Such a shell would ap- 
pear as a 1 kpc diffuse H^ emitting feature with 
brightness 



L36><10-i2iV^ 



10 ^^ergs ^ cm ^, 
(17) 

where = 1.5 arcsec is the averaged seeing. How- 
ever, this is below our detection limit (8.15 x 
10~^^ erg cm~^), and thus the H^ emission ex- 
cess expected at the bubble outer shell is unlikely 
to be detected. Approximately an order of magni- 
tude better sensitivity has to be reached to confirm 
this possible model. 

So far it appears that the remnants caused by 
the two completely different recent star-forming 
histories considered here are able to explain the 
observables in VII Zw403. In all of our models 
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however, we have also calculated the time depen- 
dent mixing of the metals freshly ejected by super- 
novae with the matter thermally evaporated from 
the outer shell. Figure 6 shows the time evolution 
of the mean inner gas metalicity (using oxygen as a 
tracer, see Silich et al. 2001) predicted for the su- 
perbubblcs in the two extreme cases able to match 
the observations of VII Zw403. Both cases show 
a superbubble interior metalicity largely different 
to that of the ISM detected in the optical regime 
(Z/sM = 0. 06.^0, see Table 1). Enormous differ- 
ences are found in the coeval case which rapidly 
acquires several times Zq. On the other hand, 
the continuous star formation, although present- 
ing high metalicity values, never surpasses Z©. 
This result could be used as a discriminator able 
to discard various possible recent histories of star 
formation in VII Zw403. However, as it is likely 
that the swept-up matter in the shell preserves 
the same metalicity as the unperturbed ISM (see 
Tenorio-Tagle, 1996), it will be difficult to distin- 
guish between the metalicities of two overlaid X- 
ray components in the instantaneous burst case. 

A close comparison of the HST observations 
and our data however has led us to realize that, al- 
though UV photons may escape their production 
centers and HII regions causing the extended dif- 
fuse ionized gas, the mechanical energy from the 
various star-forming centers has not had sufficient 
time to drive a large-scale remnant. This has not 
happened, despite the 5-10 Myr of evolution of 
the various subgroups in the galactic nucleus. It 
seems that the five star-forming centers are suffi- 
ciently distant from each other, that their mechan- 
ical energy is still currently being used to build 
up individual local shells. The largest He, shell 
around star cluster 4 with a radius ^ 250 pc is not 
yet enclosing, and not even in contact with the HII 
regions and smaller shells produced by the other 
star forming centers present in the galaxy. 

This structure thus invalidates the assumption 
often made, that a given mechanical luminosity 
derived from the H^ luminosity of a galaxy can be 
directly used in theoretical models as acting upon 
a smooth medium from the start of the calculation. 

VII Zw403 is clearly indicating that there may 
be an important leakage of UV photons out from 
the centers of star formation, causing an extended 
diffuse ionized region around a star burst. How- 
ever, for the mechanical energy of the various well 



separated centers of star formation to start work- 
ing together in the production of a large-scale rem- 
nant, a longer time scale is required. A sufficient 
period needs to allow for the multiple shock waves 
to merge and overlap with each other, clearly im- 
plying a longer time the larger the separation be- 
tween subgroups and the denser the ISM may be 
around the various star-forming centers. 

This time scale is particularly relevant when 
one tries to match the remnants produced by the 
stellar activity in a badly (spatially) resolved or 
distant galaxy, for which one can nevertheless de- 
rive a good estimate of its Hq luminosity. If for 
example one will ignore the detailed structure of 
the ionized gas in the VII Zw403 and use the me- 
chanical energy input rate derived directly from 
the total Hq luminosity, then little difference will 
be found between the coeval and the continuous 
star formation models. To fit simultaneously the 
larger Hq, luminosity (4.7 x 10^° erg s~^) and the 
X-ray luminosity derived from the PSPC data, 
both models would require very similar star cluster 
masses (1 — 3) x 10^ Mq, a similar ISM structure 
with a large cloud mass filling factor (~ 95%), and 
a similar superbubble dynamical age (5-8 Myr). 
In both cases the expanding shell would remain 
adiabatic when the remnant reaches the 1 kpc ra- 
dius. That is, if one takes the total H^ luminosity 
to derive the mechanical energy input rate and the 
bubble time evolution without knowing the ionized 
gas spatial distribution, one would predict a very 
powerful and young remnant. The delay caused by 
depositing the derived total mechanical luminos- 
ity locally into the different knots, clearly changes 
the outcome. Nevertheless, the high level of en- 
ergy deposition (Lsb ~ 2 x 10^° erg s^^) predicts 
(see Silich & Tenorio-Tagle, 2001) that the ISM of 
VII Zw403 will be blown away after a few tens of 
Myrs. 

This leads to another possible interpretation of 
the observed X-ray emission, that is coming from 
the local bubbles generated by the various stellar 
subgroups. Take for example the OB associations 
1 and 4 for which we know the total stellar mass, 
the ionizing flux, the age, and the local bubble ra- 
dius. From these we have derived both analytic 
(Chu & Mac Low, 1990), and numerical (Silich 
et al., 2001) estimates of the present bubble X- 
ray luminosity. Note that local bubble kinematic 
ages estimated from the Hq, radii and expansion 
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velocities (~ 1 Myr) arc much smaller than the 
stellar cluster ages derived from the massive star 
isochrones, which are 4-6 Myrs (Lynds et al. 

1998) . This discrepancy is often observed in the 
LMC bubbles (Oey & Smedley, 1998), and most 
probably is related to the recent bubble blowout 
from the host molecular cloud (Silich & Franco, 

1999) . We have assumed a mean stellar cluster 
age (5 Myrs) as a better indicator of the bubble 
time evolution. One can then estimate the X-ray 
emission from the local bubbles if the energy in- 
put rate and the surrounding gas density nc are 
known: 

= lQ^^ZI{T)Lff\^J'^Hf'''\ (18) 

where Z is the hot X-ray emitting gas mctalicity, 
I(t) is a dimensionless integral whose value is close 
to unity, the energy input rate Lsg, and the evo- 
lutionary time ty, measured in 10^* org s^^ and 
10^ yr units respectively. We derive the appro- 
priate star cluster mass, and its mechanical lumi- 
nosity from our SB model, and use them to esti- 
mate the surrounding gas density from the stan- 
dard (Weaver et al. 1977) bubble model: 



/267pcY 

V Rsh J 



381-7, 



(19) 



These results were compared (sec table 4) with the 
numerical calculations assuming a time-dependent 
mechanical energy input rate and time-dependent 
hot gas mctalicity. Note that the inferred X-ray lu- 
minosities show a good agreement with the obser- 
vations. The large differences between the numer- 
ical and analytical estimations of the surrounding 
gas density results from the different mechanical 
energy input rates (a constant value for analytical 
models and increasing with time rate for numerical 
calculations). 

6. Conclusion 

Here we have discussed the recent history of 
star formation and a possible nature of the dif- 
fuse X-ray emission in the nearby BCD galaxy VII 
Zw403. Two possible scenarios of star formation 
have been considered: an instantaneous burst, and 
an extended episode of star formation. 

To construct the numerical model we have pro- 
vided new narrow-band observations of VII Zw403 



centered on the Hq, line with a long exposure time 
corresponding to the threshold Hq, flux 8.15 x 
10~^^ erg cm~^ s~^. These observations reveal a 
variety of ionized filaments and incomplete shells 
superimposed on the diffuse emission that 
most certainly result from the photons leaking out 
of the main star-forming centers. The largest fea- 
ture is the 250 pc broken shell associated with stel- 
lar association 4. The total luminosity derived 
from our observations, Lh^ = 4.7 x 10^" erg s~^, is 
much larger than reported before. This has a pro- 
found impact on the predicted properties of the 
starburst blown superbubble. 

The numerical models based on the HST Ra 
data require either an instantaneous burst of star 
formation with a total mass of 5 x 10^ Mq , or a star 
formation episode with a constant SFR = 4 x 10~^ 
Mq yr~^ lasting 35 Myr. The models however re- 
quire radically different structures of the galactic 
ISM and imply very different properties of the re- 
sulting remnant. 

The best coeval model assumes most of the ISM 
to be locked up within high density clouds, and 
only ~ 5% of the observed neutral hydrogen mass 
is in the smooth component. The hydrodynami- 
cal calculations also predict the outer shell to be 
adiabatic after reaching a 1 kpc radius, and to 
contribute 50 - 80% of the observed diffuse X-ray 
emission. The bubble evolutionary time is esti- 
mated to be Tdyn ~ 7 Myr when its expansion 
speed is w 200 km s~^. 

The best continuous star formation model re- 
quires of much higher density in the smooth ISM 
component, with only ~ 30% of the HI mass con- 
centrated in dense clouds. This leads to a much 
smaller bubble expansion velocity (Vexp ~ 35 km 
s~^), larger evolutionary time (rdyn ~ 35 Myr), 
and a rapid cooling within the outer shell. That 
is, this model predicts a low brightness HI shell 
surrounding the diffuse 1 kpc X-ray region. This 
shell may also show up in Ha if exposed to the 
UV flux from the central cluster. The inner gas 
metalicities are also predicted to be very different 
in these two cases. 

The numerical calculations based on the high 
energy input rate derived from our observations 
require an instantaneous burst or a short episode 
of star formation with SFR ~ 0.1 Mq yr~^ lasting 
less than 10 Myr with similar total stellar cluster 
masses (1-3) x 10^ Mq and most of the ISM 
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(~ 95%) locked up within high density clouds. 
The comparison of the energy input rate derived 
from our Ha data with the theoretical limits, im- 
plies that the entire ISM and metals produced by 
the current episode of star formation are going to 
be ejected from the galaxy after the coherent su- 
per bubble is formed. 

It appears that the five main star-forming knots 
are sufficiently distant to form a coherent shell in 
a short time scale, while keeping their energies 
blocked within local, spatially separated bubbles. 
This provides a time delay that must be considered 
when developing a numerical model for the coher- 
ent superbubble driven by a number of young stel- 
lar clusters. Numerical calculations show that the 
X-ray luminosities from young local bubbles are 
in a good agreement with the ROSAT PSPC data. 
This agreement indicates that the observed diffuse 
component of the X-ray emission may be related 
to the small centrally concentrated bubbles, rather 
than to the coherent 1 kpc structure. Further ob- 
servations with the XMM-NEWTON observatory 
is expected to be able to recover the real nature of 
the diffuse X-ray emission and the recent history 
of star formation in this galaxy. 
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Table 1: VII Zw403 observational properties 



Parameter 




Keierence 


Comments 


Distance 


4.0 Mpc 


-1 
1 




ine galaxy total mass 


^ X ID M0 


2 




Lhe galaxy IbM rticiss 


(4 — 7) X 10 M© 


1, 2 


from HI observations 


The radius of tlie HI halo 


~ 1.9 kpc 




not well known, see text 


-L lie iOiVi gcLb Velocity 








dispersion 


- 30 km s"^ 


2 




The ISM gas metaUcity 


(0.05 - O.O6)Z0 


3 




The total Ha flux 


(195 ± 4) X IQ-i^ ergs cm'^ s"^ 


our data 


integrated out to 




1.8 X 10^^ ergs s"^ 
4.7 X 10''° ergs s"^ 




the limiting isophote 


The total Ha luminosity 


1 

our data 




Ha shell radii 








- Large 


~ 250 pc 


our data 


around association 4 


- Small 


~ 80 pc 


1 


around association 1 


Diffuse X-rajr emission: 


(1.9 - 2.3) X 10^* ergs s"^ 


4, 5 


ROSAT PSPC data 


- Unresolved core 


-85% 




not confirmed by 


- Extended diffuse emission 


~15% 




ROSAT HRI data 


X-ray flux from the point 








source 


9 X 10"^^ ergs cm"^ s"^ 


6 


ROSAT HRI data 



References: 1) Lynds et al., 1998 2) Thuan & Martin, 1981 3) Schulte-Ladbeck et al., 1999 4) Papaderos et al., 1994 5) Fourniol, 1997 
6) Lira et al., 2000 



Table 2: A 1 kpc bubble parameters for instantaneous burst models 



Model 


MsB 


fc 


'^dyn 


Nuv 


Ljfa 




Shell state 




Mo 


% 


Myr 


photons s~^ 


er 


g S 


erg s~^ 




1 


2 


3 


4 


5 


6 


7 


8 


15.0.10 


lO'' 


90 


12.5 


7.4x10'*'' 


1.0 


X 10'" 


1.2 xlO-" 


R 


15.0.40 


10'' 


60 


21.8 


4.3x10"^ 


5.8 


xlO^" 


2.2 xlO^^ 


R 


15.0.70 


10'' 


30 


29.7 


8.7x10'"' 


1.2 


xlO^'^ 


2.8 xlO^^ 


R 


15.7.5 


5 X lO'"' 


95 


6.7 


9.5x10-'''" 


1.3 


xlO-^» 


4.2 xlO-^^ 


A 


15.7.10 


5 X lO'"' 


90 


7.9 


4.3x10™ 


5.8 


X 10-^** 


1.6 xlO-^*^ 


A 


15.7.40 


5 X lO'"' 


60 


12.0 


4.5x10'*^ 


6.1 


X 10'^^ 


1.8 xlO^* 


R 


15.7.70 


5 X lO'"' 


30 


15.5 


1.3x10'*^ 


1.7 


xlO^'^ 


2.3 xlO^* 


R 


16.0.5 


lO'' 


95 


6.0 


3.5x10^^ 


4.8 


xlO*^ 


3.5 xlO^^ 


A 


16.0.10 


10'' 


90 


6.8 


2.0x10"^ 


2.7 


xlO-^'' 


1.5 xlO^* 


A 


16.0.40 


lO** 


60 


9.8 


2.7x10-'''" 


3.7 


X 10'^** 


4.5 xlO^** 


R 


16.0.70 


lO** 


30 


12.0 


9.1xl0''« 


1.1 


xlO^** 


5.5x10-^** 


R 


17.0.10 


10^ 


90 


4.7 


1.5x10"^ 


2.0 


xlO^i 


3.3x10=^® 


A 


17.0.40 


10^ 


60 


5.8 


4.4x10^^ 


5.9 


xlO« 


1.7x10^^ 


A 


17.0.70 


10^ 


30 


6.3 


3.0x10^^ 


4.1 


xlO« 


8.1x10^^ 


A 
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Table 3: A 1 kpc bubble parameters for continuous star formation models 



Model 


MsB 


fc 




Nuv 




La; 


Shell state 




Mo 


% 


Myr 


photons 


erg s^^ 


erg s^''" 




1 


2 


3 


4 


5 


6 


7 


8 


C5.10 


2.1x10* 


90 


22.3 


1.3x10*' 


1.6x10^'' 


1.3x10^" 


R 


C5.40 


2.1x10* 


60 


50.5 


1.5 xlO*'^ 


2.1x10'^'^ 


1.1x10^® 


R 


ClO.lO 


4.1x10* 


90 


20.0 


1.5x10** 


2.0x10'^'* 


3.6x10^'^ 


R 


C10_40 


4.1x10* 


60 


34.5 


4.0x10*"^ 


5.4x10^* 


6.6xl03« 


R 


C10.70 


4.1x10* 


30 


59.0 


1.7x10*^ 


2.3x10^^ 


no 


R 


C20.10 


8.2x10* 


90 


19.5 


1.3x10^^^ 


1.8x10^'' 


5.8 X 10^*^ 


R 


C20_40 


8.2x10* 


60 


29.4 


1.7x10*** 


2.3x10'^*' 


1.6x10*^ 


R 


C20.70 


8.2x10* 


30 


37.5 


1.6x10*^ 


2.2x10'^" 


4.3x10^^ 


R 


C40.10 


1.6x10'' 


90 


19.5 


1.3x10'''^ 


1.8x10^^ 


6.3x10-*^ 


R 


C40_40 


1.6x10^ 


60 


29.0 


1.3x10^^ 


1.8x10*^ 


1.9x10^^ 


R 


C40.70 


1.6x10^ 


30 


35.1 


1.3x10^^^ 


1.8x10^'-' 


3.2x10^^ 


R 



Table 4: Parameters from individual bubbles 





Association 1 


Association 4 


UV flux 


5 X 10'^" s-i 


7 X lO^'^ s-i 


Star cluster age 


5 Myr 


5 Myr 


Star cluster mass 


4.7 X 10"* M0 


6.6 X 10'' M0 


Shell radius 


79 pc 


250 pc 




Analytic model 


Analytic model 


ISM number density 


1.1 X 10-* cm-^ 


4.9 cm~3 


Bubble X-ray luminosity 


3.5 X 1038 Z I(t) erg s"^ 


3.4 X 10^7 Z I(t) erg s'^ 




Numerical model 


Numerical model 


ISM number density 


40 cm~3 


0.2 cm-3 


Bubble X-ray luminosity 


1.8 X 10^8 erg s'^ 


2.1 X lO^^^ erg s"^ 


Hot gas mean metalicity 


Zo = 1.8^0 


Zo = O.GZq 
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Fig. 1. — The intrinsic properties of instantaneous and extended starbursts. The panels present: the energy 
input rates (a), the ionizing photon flux (b), and the total mass, and the fraction of this in oxygen, ejected 
by instantaneous and extended bursts of star formation (c), all as a function of time. In all the panels solid 
lines represent the instantaneous {tsb = 1 Myr) burst, dotted lines indicate the results for the tsb = 20 
Myr and dashed lines those for r^s = 50 Myr. 
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Fig. 2. — The narrow-band Ha image of VII Zw403. The left panel shows isocontours map of the continuum 
substracted Ha image. The right panel presents grey scale map of Hq emission in logarithmic scale. The 
lowest Hq, isocontour level corresponds to the threshold value 8.15 x 10"^'' erg s~^ cm~^. Isocontours are 
equispaced in logaritmic scale. 
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Fig. 3. — The gas density distribution for models assuming 90% of the ISM mass to be stored in clouds. 
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Fig. 4. — Inferred parameters for VII Zw403. Prom the observed Hq, luminosity one can infer the SFR and 
thus the ionizing photon flux. If all those quantities remain constant for a tsb as shown in panel (a) for the 
central region Hq, luminosity, one can then infer the run of the mechanical energy input rate (b). 
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Fig. 5. — Comparison of the X-ray luminosity (diffuse component) for models with tsb = 40 Myr and 
different fractions of the interstellar gas in the cloudy component (Mgi). 
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Fig. 6. — The metalicity of the superbubble as a function of time for two best models related with the HST 
data. Oxygen is used as a tracer, arrows indicate metalicities at the time super bubbles reach the 1 kpc 
radius. 
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